The effect of limb kinematics on the speed of a legged robot on granular media 
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Achieving effective locomotion on diverse ter- 
restrial substrates can require subtle changes of 
limb kinematics. Biologically inspired legged 
robots (physical models of organisms) have shown 
impressive mobility on hard ground but suffer 
performance loss on unconsolidated granular ma- 
terials like sand. Because comprehensive limb- 
ground interaction models are lacking, optimal 
gaits on complex yielding terrain have been de- 
termined empirically. To develop predictive mod- 
els for legged devices and to provide hypothe- 
ses for biological locomotors, we systematically 
study the performance of SandBot, a small legged 
robot, on granular media as a function of gait 
parameters. High performance occurs only in a 
small region of parameter space. A previously in- 
troduced kinematic model of the robot combined 
with a new anisotropic granular penetration force 
law predicts the speed. Performance on granular 
media is maximized when gait parameters mini- 
mize body acceleration and limb interference, and 
utilize solidification features of granular media. 



INTRODUCTION 

To move effectively over a wide range of terrestrial 
terrain requires generation of propulsive forces through 
appropriate muscle function and limb kinematics [TJ 
2 . Most biological locomotion studies have focused 
on steady rhythmic locomotion on hard, flat, non-slip 
ground. On these surfaces kinematic (gait) parameters 
like limb frequency, stride length, stance and swing dura- 
tions, and duty factor can change as organisms walk, run, 
hop and gallop lj. There have been fewer biological stud- 
ies of gait parameter modulation on non-rigid and non- 
flat ground, although it is clear that gait parameters are 
modulated as the substrate changes during challenges like 
climbing OE], running on elastic/damped substrates [5], 
transitioning from running to swimming [6] , and running 
on different preparations of granular media [7J. Even 
subtle kinematic changes in gait can lead to major dif- 
ferences in limb function [5]. A major challenge is to de- 
velop models of limb interaction with complex substrates 
and to develop hypotheses for how organisms vary gait 
parameters in response to substrate changes. 



The RHex class of model locomotors (robots) has 
proved useful to test hypotheses of limb use in biolog- 
ical organisms on hard ground [9] and recently on more 
complex ground with few footholds (TU] or the ability to 
flow [TT]. These hexapedal devices model the dynami- 
cally stable locomotion of a cockroach and were the first 
legged machines to achieve autonomous locomotion at 
speeds exceeding one body length/s. In these devices, 
complexity in limb motion is pared down to a few biologi- 
cally relevant parameters controlling intra-cycle "stance" 
and "swing" phases of 1-dof rotating limbs (referred to 
as "gait" parameters hereafter; see detailed description in 
Methods and Results). When these gait parameters are 
appropriately adjusted, RHex shows performance compa- 
rable in speed and stability to organisms on a diversity 
of terrain [T^]. However, because of the scarcity of exist- 
ing models of limb interaction with complex substrates, 
adjustment of the gait parameters is typically done em- 
pirically USE]. 

Sand, a granular medium |15j . is of particular interest 
for studies examining the effects of limb kinematics on lo- 
comotor performance on yielding terrain. In a previous 
study [H] we found that minor changes in the limb kine- 
matics of a small RHex-class robot, SandBot (Fig. la), 
produced major changes in its locomotor mode and per- 
formance (speed) on a granular medium, poppy seeds 
(see Section 2.2). This sensitivity occurs, in part, be- 
cause forced granular media remain solid below the yield 
stress, but can flow like a fluid when the yield stress is 
exceeded [IB]- We tested SandBot on granular media of 
different yielding properties (set by granular volume frac- 
tion; see Section 2.1) at various limb frequencies but with 
the other gait parameters fixed. While there is no fun- 
damental theory at the level of fluid mechanics that ac- 
counts for the physics of the solid-fluid transition of gran- 
ular media or the dynamics of the fluidized regime, em- 
pirical models of granular penetration force have proved 
useful to predict SandBot 's speed [11]. SandBot 's propul- 
sion is determined by factors that control this transition 
during limb-ground interaction (limb penetration depth, 
limb speed, body mass, grain friction, volume fraction, 
etc.). Using a simplified equation describing the granu- 
lar penetration force, we developed a kinematic model to 
explain the locomotion of SandBot (see Section 2.2). 

In this study, we advance our understanding of the 
effects of limb kinematics on locomotor performance by 
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testing SandBot with varying gait parameters on sand of 
fixed yield strength and at fixed limb frequency. We find 
that robot speed depends sensitively on limb kinemat- 
ics; while the original model qualitatively captures this 
sensitivity, the penetration force used in the model and 
other assumptions need to be modified to explain some 
important features. Our study not only reveals the spe- 
cific optimal kinematics for SandBot on granular media, 
but also advances our understanding of how in general to 
achieve effective legged locomotion on complex terrain. 



BACKGROUND AND REVIEW OF PREVIOUS 
STUDY 



cycle- averaged limb frequency u>. We observed a sensi- 
tive dependence of v x on both <j> and u, and developed a 
kinematic model which explained this dependence and 
revealed two distinct locomotor modes determined by 
whether the granular media solidifies during limb-ground 
interaction. 

Our kinematic model describes the limb-ground in- 
teraction of SandBot by considering the motion of just 
a single limb (Fig. lb). SandBot has six approxi- 
mately c-shaped limbs (c-legs) divided into two alternat- 
ing tripods. C-legs in the same tripod rotate in synchrony 
and each c-leg rotates about a horizontal axis normal to 
the robot body. We simplify the multi-leg ground inter- 



To understand the effect of limb kinematics addressed 
here, we first summarize the mechanism of SandBot lo- 
comotion on granular media (called rotary walking) dis- 
covered in our previous study In this section, we 
discuss the physics of granular media that controls the 
limb penetration depth (which governs locomotion per- 
formance) and then review our previous experiments and 
kinematic model. 



Physics of Limb-Granular Media Interaction 

The physics that controls locomotor performance is 
the relative magnitude of the penetration resistance force 
(originating in the granular media) and the sum of the 
external forces (weight, inertial forces). When these bal- 
ance, the granular media solidifies, allowing the robot to 
be supported at a fixed limb penetration depth. 

The previous SandBot study [11] revealed that as the 
limb (or any simple intruder) vertically penetrates into 
the medium, the penetration force scales with z, the 
depth of the intruder below the surface [T7], as F p (z) — 
k(4>)z, where </> is the volume fraction, the ratio of the 
solid volume of the granular media to the volume that it 
occupies (for natural dry sand, 0.55 < <j) < 0.64). The 
constant, k((f>), characterizes the penetration resistance 
and increases with (f>. In this paper we keep <fi fixed at ap- 
proximately the critical packing state [TH1 HH HH] (which 
is close to the as-poured volume fraction) where granular 
media neither globally dilate nor compact in response to 
shear (see Section 3). 



Review of Previous Observations and Model 

In the previous study of SandBot [11], we fixed intra- 
cycle limb kinematics (by using gait parameters that 
produce consistent motion on granular media, see Sec- 
tion 3) and measured SandBot's average speed v x on 
poppy seeds as a function of volume fraction and the 
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FIG. 1. Mechanism of SandBot locomotion on granular me- 
dia, (a) SandBot, a six-legged insect inspired robot, moves 
with an alternating tripod gait. The three arrows indicate 
the limbs of one tripod, (b) Schematic of single-leg represen- 
tation of SandBot with mass m — 1/3 SandBot's total mass. 
With the body contacting the surface, the motor axle is height 
h — 2.5 cm above the ground. The c-leg is approximately a 
circular arc (radius R = 3.55 cm, arc span 225 degree). Leg 
angle 9 is measured clockwise about the axle and between 
the downward vertical and a diameter through the axle. Leg 
depth z = 27? cos 9 — h. (c) Magnitude of penetration force 
F p (blue curve) relative to force required for upward motion, 
m(g + a), (red lines) determines the locomotor mode. The 
force required for quasi-static movement (mg) is shown for 
reference. When F p and m(g + a) intersect rotary walking 
occurs. When F p and m(g + a) do not intersect (dashed red 
curve and above), the robot swims, (d) Schematic of rotary 
walking. The granular material flows in the intervals [on, j3i\ 
(red arrow) and [/3/, q/] (blue arrow) where F p < m(g+a) and 
the c-leg rotates about the axle (red and blue circles). The 
material is a solid in the interval [/3i, /3/] (gray sector; line 
with arrows in (c)) where F p exceeds m(g + a) and the c-leg 
rotates about its center (green circle and arrow), lifting and 
propelling body forward by step length s — i?(sin/3/ — sin/3j). 
The c-leg is above the ground in the interval [ctf, on + 2ir]. 
Note that [at , a/] in (c) is symmetric to vertical (9 = 0) as a 
result of assuming the force is isotropic (see Section 4.2). 
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action of each tripod to that of a single limb carrying 
1/3 the total body mass 3m (2.3 kg), as body weight is 
approximately uniformly distributed between each c-leg. 
We also considered a c-leg as a simple intruder ignoring 
its more complicated geometry, i.e. F p (z) = kz. The 
previous study showed that the simple intruder ap- 
proximation gave approximately the same results as a 
more realistic treatment in which penetration force was 
integrated over the submerged leading surface of a c-leg. 
In this study we use the simple intruder approximation. 

We make the approximation that SandBot's body is in 
stationary contact with the surface at the onset forward 
motion in each cycle [20], and define the c-leg's angular- 
position, 9, as the clockwise angular displacement from 
the configuration where the center of curvature of the c- 
leg is directly beneath the axle, see Fig. lb. During a 
full rotation, as 8 changes from — ir to n, the c-leg ini- 
tially contacts the ground at 9 = a, and loses ground 
contact at 6 = a t . Because leg depth can be approxi- 
mated as z — 2R cos 9 — h when the body is in contact 
with the surface [UJ , penetration force can be written as 
F p (9) — 2Rk [cos 9 — ^] (blue curve in Fig. lc), where 
R = 3.55 cm and h = 2.5 cm are the radius of the c- 
leg and the hip height (i.e. distance from c-leg axle to 
underside of body) respectively. 

Of prime importance in determining SandBot's per- 
formance is the magnitude of the penetration force F p (9) 
relative to the sum of the forces required to the sup- 
port the body weight and accelerate the body upward 
m(g + a) (red curve in Fig. lc), where g is the accel- 
eration due to gravity and a the acceleration [TT. The 
relevant acceleration is given by the jump in robot speed 
when the granular media solidifies, Rlu, divided by the 
characteristic response time of the c-leg interacting with 
the granular media, At((f>), i.e. a — Rui/At. Two distinct 
locomotor modes are possible depending on whether or 
not F p (6 = 0) > m(g + a): 

1. Rotary walking - movement with solidification (see 
Fig.lb-d): As the c-leg rotates into the ground after 
initial leg-ground contact at 9 = cti, the penetra- 
tion force increases with increasing depth. In the 
rotary walking regime the material beneath the c- 
leg solidifies and leg penetration stops at an angle 
9 = (3i when F p (/3i) = m(g + a), see Fig. lc,d. Since 
the frictional force between the c-leg and granu- 
lar material is insufficient for the leg to roll, the 
c-leg instead rotates about its center of curvature 
(green circle and arrow) lifting and advancing the 
robot in the process. Rotary walking continues 
until 8 = (3f, beyond which the c-leg again pene- 
trates through the material since F p (9) < m(g + a) 
and the body is again in contact with the ground 
(blue circle and arrow). Rotary walking thus oc- 
curs over a finite range of leg angle < 9 < (3f or 
[/?i,/3/] (horizontal arrow in Fig. lc and gray sec- 



tor in Fig. Id) where fii and /?/ are determined 
by F p {pij) = 2Rk(cos/3ij - ^) = m{g + a). For 
a given [/3i,/3f], Fig. Id shows that the robot ad- 
vances a distance s = R(sm(3f — sin^), where we 
call s the step size. During one complete gait cy- 
cle of period T, each alternating tripod advances 
the robot by s, giving an average robot speed of 
v x = 2s/T = suj/tt. 

2. Swimming - movement without solidification: 
When F p (0) < m(g+a) (Fig. lc, dashed red curve), 
the granular material beneath the penetrating c-leg 
never solidifies and rotary walking does not occur, 
i.e. /3j = j3f — 0. Instead, the limb constantly slips 
through the surrounding fluidized granular mate- 
rial, similar to a swimmer's arm in water, and the 
robot advances slowly (v x < 1 cm/s). In this 
regime forward motion occurs when the frictional 
and inertial (drag) forces generated by the c-legs 
exceeds the frictional force between the robot body 
and the surface. 

The two constants characterizing the interaction with 
the granular medium, k and At, together with limb fre- 
quency w, determine the relative magnitudes of F p and 
m(g + a) and consequently control which locomotor mode 
the robot operates in. Reducing k (by decreasing (j>) 
and/or increasing lu reduces the rotary walking range; 
in other words, the less compact the granular material is 
and/or the faster the limbs rotate, the deeper the c-legs 
have to penetrate before the granular material solidifies 
and rotary walking begins, and the more susceptible the 
robot is to entering the slow swimming mode. This sim- 
ple kinematic model captures the observed sensitive de- 
pendence of v x on <p and u>, with k(<p) and At(<j)) as two 
fitting parameters. 

In summary, our previous study of SandBot [TT] 
showed that to locomote effectively on granular media, 
limbs kinematics that access the solid phase of granular 
media should be employed. 

METHODS AND RESULTS 

The limb kinematics of each tripod during one cy- 
cle are parameterized by three "gait parameters", see 
Fig. 2(a,b). The kinematics of both tripods are periodic 
(with period T) and offset by half a period T/2 but are 
otherwise identical. For the conditions in this and previ- 
ous experiments, a motor controller in the robot ensures 
that the target kinematics are achieved. Limb kinemat- 
ics consist of a "swing" phase (orange), which is typically 
faster, and a "stance" phase (green), which is typically 
slower, with respective frequencies uif and ui s . 

During hard ground locomotion in the RHex-class of 
Robots (and for animal locomotion in general), "swing" 
and "stance" phases typically correspond to off-ground 
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FIG. 2. SandBot's intra-cycle limb kinematics and affects 
its speed on granular media, (a) Each leg rotation is com- 
posed of a fast phase (orange) and a slow phase (green). 
8 S and 80 define the angular extent and center of the slow 
phase respectively, (b) Leg angle 8 as a function of time 
during one cycle (normalized to T). 8(t) of the other tri- 
pod is shifted by T/2 but otherwise identical. d c is the 
duty cycle of the slow phase, i.e. fraction of the period 
spent in the slow phase, (c) Instantaneous speed of Sand- 
Bot on granular media with hard ground clock signal (HGK: 
{8 B ,do,d c } = {0.85,0.13,0.56}; red) and soft ground clock 
signal (SGK: {8 s ,8 ,d c } = {1.10,-0.50,0.45}; blue). With 
HGK SandBot moves slowly (v x w 2 cm/s) on granular me- 
dia, but with SGK (red), it advances rapidly (v x « 8 cm/s). 



and ground-contact phases, respectively. But because 
during locomotion on granular media this correspondence 
is not necessarily true, we simply call them fast and slow 
phases. In practice the fast and slow phases are implicitly 
defined by the triplet {9 s ,9q, d c } where 9 S is the angular 
extent of the slow phase, 9$ is the angular location of 
the center of the slow phase, and d c is the duty cycle 
of the slow phase (the fraction of the period in the slow 
phase). Specifying the cycle averaged limb frequency ui 
fully determines the motion of the limbs in the robot 
frame. By definition, ui s — ^j-, L)f = ) ' an< ^ w = 

-jr. Typically, gait parameters are set so that ui s < ui < 
LOf, but the reverse is possible when 8 S becomes large 
enough and/or d c small enough. 

In the first tests of SandBot on granular media 
(Fig. 2c), we found that kinematics tuned for rapid stable 
bouncing motion on hard ground (HGK: {9 s ,0o,d c } — 
{0.85, 0.13, 0.56}) produced little motion on granular me- 
dia (red curve). Empirical adjustment to soft ground 
kinematics (SGK: {6 s ,6 ,d c } = {1.10,-0.50,0.45}) re- 
stored effective (walking) locomotion on granular media 
(blue curve). 



In the previous study [TT], we used SGK to test 
v x (<j),ui). Now armed with the understanding of how 
SandBot moves on granular media gained from this work, 
we set out to determine the effects of limb kinematics in 
detail. We set <p = 0.605 and u = 8 rad/s, and measure 
SandBot's average speed on granular media as we sys- 
tematically vary gait parameters, i.e. v x — v x (9 s ,9o,d c ). 
We pick u> = 8 rad/s because at this intermediate fre- 
quency SandBot displays both rotary walking and swim- 
ming as the clock signal is varied. We pick = 0.605 to 
remove the effect of local volume fraction change which 
causes a premature transition from rotary walking to 
swimming [22j and adds to the complexity of the prob- 
lem. 

We first test the effect of the extent and location of 
the slow phase for fixed d c = 0.5, measuring speed v x = 
v x (9 s ,9q). We vary the parameters between < 9 S < 2 
and —2<9o< 2, which are the limits set by the robot's 
controller. We choose d c = 0.5 because it is close to the 
d c values of both HGK and SGK. This gave us an easy 




FIG. 3. Average speed v x of SandBot on granular media 
(4> — 0.605) as a function of (8 3 , 8q) for d c = 0.5 and cu = 
8 rad/s. (a) The experimental data has a localized region 
of high speeds with peak v x ~ 9 cm/s near {8 3 ,8o} = {1.5, 
-0.5}.. Circles show that v x for SGK (red) and HGK (blue) 
matches data for d c — 0.5 despite the formers slightly different 
d c values, see Section 4.3. Inset: original data from which 
main figure is interpolated, (b) Predicted v x (8 3 , 80) from the 
kinematic model with F p = kz captures the single peak but 
predicts a lower speed for SGK than for HGK contrary to 
observation and fails to account for the observed peak in speed 
at 8 = -0.5. 
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way to project HGK (d c = 0.56) and SGK (d c = 0.45) 
onto the v x = v x (9 s ,9q) plot (d c = 0.5), assuming that 
a small change of d c near d c = 0.5 does not affect speed 
significantly (see Fig. 6a and Section 4.3 which support 
this assumption)). 

Measurements of v x = v x (9 s , 9q) (Fig. 3a) show a sin- 
gle sharp peak in speed near {9 S , 9q} = {1-5, —0.5}. High 
speeds only occur within a small island of — 1 < 9q < 
and 9 S > 0.5 surrounding the peak; lower speeds fill the 
remainder of the space. The drop in speed is rapid as 
#o is varied away from the peak, and is less so when 
S is varied away from the peak; this is also evident in 
cross sections through the peak (blue circles in Figs. 3a 
and Fig. 6a, respectively). Ignoring the effect of d c , the 
SGK parameters (blue dot) lie close to the peak while 
the HGK parameters (red dot) are in the low speed re- 
gion. The optimal gait parameters which we found for 
SandBot locomotion on poppy seeds at (f> = 0.605 and 
w = 8 rad/s are: {6 s ,6 ,d c } = {1.5,-0.5,0.55}. These 
gait parameters generate about 20% higher speed than 
the previously used SGK parameters. 

Variation of the duty cycle at fixed {9 S ,9 } = 
{1.5,-0.5} also has a substantial influence on speed. 
Data (blue circles in Fig. 6b) show a well defined peak 
at d c « 0.5. Speed drops off relatively slowly for d c > 
0.55, and more quickly to small (swimming) speeds for 
d c < 0.5. 

DISCUSSION 

Application of Model to Slow Phase Extent and 
Location Variation 

To apply our kinematic model to SandBot locomotion 
with varied limb kinematics, we must consider the effects 
of variable limb kinematics during limb-ground interac- 
tion. Depending on the gait parameters during ground 
contact, the limbs could be rotating in the fast phase, in 
the slow phase, or in a combination of both. In our pre- 
vious study, the kinematic model ignored limb frequency 
variability during ground contact and only considered the 
robot limb rotating at the constant cycle averaged limb 
frequency uj. 

However, as limb kinematics change, the variability of 
limb frequency in ground contact needs to be taken into 
account. For our test of v x = v x (9 s , 9 ) within < 9 S < 2 
and — 2 < 9 < 2 and at d c = 0.5, uif >> lu s so that only 
the slow phase can possibly achieve rotary walking, as 
fast limb rotation results in swimming. In this case oj b 
(instead of w) controls acceleration a and thus determines 
the rotary walking range |23) . 

For fixed d c , varying [9 s ,9q] changes the extent and 
location of the slow phase; the angular extremes of the 
slow phase are 9ij — 9 Q ± ^f, see Fig. 5a. Varying 9 S 
also changes oj s which controls the rotary walking range. 



Therefore varying [0 s ,#o] affects where the slow phase 
overlaps with the rotary walking range. The step length 
s is given by s = R(smtjjf—smtpi) where ipi = max(ft, 9i) 
and tjjf = min(/3j, 9f) if there is overlap or s = if there 
is no overlap. The larger the overlap, the further the 
robot moves forward in a cycle. 

As shown in Section 2.3, the rotary walking range 
[/3i,f3f] is given by solving the equation F p {9) = 
2Rk(cos/3ij- ^) = m(g + a), with a given by a = 
We can evaluate how [9i,9f] overlaps with to de- 

termine s, and calculate the robot speed using v x — |f — 
*jjr. For fixed uj, speed v x scales with step length s. 

Figure 3b shows the model prediction of v x using fit- 
ting parameters k = 210 N/m, At = 0.37 s. Compar- 
ing prediction with observation (Fig. 3a), the model cap- 
tures the peak and predicts similar magnitudes of speeds. 
However the predicted peak is symmetric about 9q = 
while the observed peak is symmetric about 9 = —0.5. 

Anisotropic Penetration Force Law 

If the penetration force of the granular material in- 
creased like kz as assumed in the model, we would ex- 
pect = as this value would give the largest overlap 
between the slow phase and the rotary walking range 
as determined by the material strength (see scheme in 
Fig. 5a) . To investigate why the robot performs best with 
#o = —0.5 we attached a c-leg to a force/torque sensor 
and measured the grain resistance as the c-leg was ro- 
tated through the granular media at uj = 0.35 s _1 (the 
horizontal rotation axis of the c-leg was positioned the 
same distance h = 2.5 cm above the grain surface as 
when it is mounted on the robot). Figure 4b shows a 
clear asymmetry in the penetration (vertical) force with 
the measured peak force occurring near fto = —0.75; pen- 
etration force during rotation peaks before the intruder 
reaches the maximum depth. We confirmed that the 
measured anisotropy in the penetration force is intrin- 
sic to our granular medium and is not an artifact of the 
particular shape of the c-leg by additionally rotating a 
rectangular bar and a sphere into granular media at the 
same hip height: both objects exhibited a peak force at 
A, = -0.75. 

We speculate that the asymmetry in penetration force 
during rotation into granular media is a result of the 
changing limb orientation during rotational intrusion. 
For vertical penetration (which we considered in the 
model in the previous study), the intruder is constantly 
pushing down on the granular material. The grain con- 
tact network generated in granular material in response 
to intrusion |24j forms a downward pointing cone which 
generates a force symmetric to the vertical (9 — 0). In 
rotational intrusion, however, the direction of intrusion 
is constantly changing; the direction of the force cone 
should change as well and correlate with the instanta- 
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FIG. 4. Asymmetry of v x with respect to 9o = is due to 
anisotropic penetration force during limb rotation into gran- 
ular media, (a) For all 9 S at d c = 0.5, v x (9,9o) (Fig. 2a) is 
maximal (dashed vertical black line) at 9q = —0.5 (9 S = 1.5 
shown). Inset: peak location 9o — —0.5 does not change for 
d c — 0.8. (b) Vertical penetration force F p (solid blue curve) 
during c-leg rotation into poppy seeds reaches maximum at 
/3q = —0.75 (dashed black line) and is asymmetric to 9o = 0. 
Inset: force measurement schematic. Solid blue curve in (a) is 
prediction from the model with anisotropic penetration force. 



neous direction of intrusion. 

We hypothesize that the force during rotational intru- 
sion is maximal at /3q = —0.75 because for larger angles 
part of the cone reaches the surface and/or terminates on 
the horizontal walls of the container and can no longer 
support the entire grain contact network, thus reducing 
the maximal yield force. We also note that the angle at 
which maximal force is developed is close to the angle of 
repose 0.52 that we measure for the poppy seeds. This 
angle is the same as the internal slip angle in cohesionless 
granular material |16) which plays an important role in 
the formation of the grain contact network, supporting 
the plausibility of our speculation. 

To account for the measured angular offset in peak 
force from vertical (Figs. 2c and 3a), we modify the orig- 



inal penetration force law in our model to 

F p {9) = 2Rk' {cos [b (0 - fa)] + 1} 

for F > 0, where (3q — —0.75, and k! and b are new 
fit parameters. Following the same procedure described 
in Section 4.1, we find the robot speed by calculating 
[(3i,flf], the overlap between and [/3j,/3j], and the 

step size. 

Figure 5c shows v x predicted by a fit to the model using 
the anisotropic penetration force law (fitting parameters 
k' = 65 N/m, At = 0.4 s, and b = 0.8). Besides cap- 
turing the peak behavior of measured speed, the model 
also captures the shift in peak location to 6q = —0.5 (i.e. 
asymmetry to 8q = 0). For fixed 6 S , speed is maximal 
when the center of the slow phase corresponds with the 
center of the rotary walking range (Fig. 5a) . If 9o is differ- 
ent from /?o = —0.75, the overlap of the slow phase and 




FIG. 5. Overlap of the slow phase and the rotary walk- 
ing interval predicted by the anisotropic penetration force 
law better predicts v m (9 s ,9o). (a) Overlap of the slow phase 
[9i,9f] and the rotary walking range [/3i,/3f] determines step 
length s and thus speed v x . For the configuration shown, 
s — i?(sin/3/ — sin 9i). 9a and flo are centers of the slow 
phase and the rotary walking range, respectively, (b) In 
SGK the slow phase (centered at 9o = —0.50) overlaps 
nearly completely with the rotary walking range (centered at 
/?o = —0.75), explaining its high speed as compared to HGK 
(centered at 9o = 0.13) which has little overlap, (c) Pre- 
diction of v x (9 s ,9o) from the the model with the anisotropic 
penetration force for d c = 0.5 captures the asymmetry of v x 
with respect to 9o — and predicts higher speed for SGK 
than for HGK. Fitting parameters: k' = 65 N/m, At = 0.4 s 
and b = 0.8. 
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the rotary walking range decreases, which reduces step 
length and thus speed. In accord with observation, SGK 
(red dot) lies near the peak while HGK (blue dot) lies in 
a region of low speeds. Figure 5a, b shows that SGK has 
higher speed than HGK because the overlap between the 
slow phase and the rotary walking range (gray sector in 
Fig. 5a) is significantly larger. 

At fixed 6*o = /?o> increasing the extent of the slow 
phase (increasing 9 S ) from zero initially increases speed 
as the extent of the slow phase increases within the ro- 
tary walking envelope (see Figs. 2a and 6a). However, 
to s increases with 9 S which increases the acceleration and 
reduces the rotary walking range. For sufficient extent 
(near 9 S = 1.5 for the data shown in Figs. 2a and 6a) the 
slow phase contains the rotary walking range and step 
length is determined by the latter. Further increase in 
9 S reduce the rotary walking range. Rotary walking is 
not possible for 9 S > B^kM ^ - A as the material is 
never strong enough to both support and accelerate the 
robot. In Fig. 6a the experimental speed is noticeably 
lower than the model prediction at the largest 9 S = 2. 
As we discuss below in regards to variation in d c , this re- 
duction is a apparently the result of tripod overlap (both 
tripods simultaneously in ground contact) which occurs 
for a greater portion of the slow phase for larger 9 S . 



Effect of Duty Cycle 

While the model prediction of v x (d c ) (blue curve in 
Fig. 6b) matches the magnitudes of speeds at intermedi- 
ate d c « 0.5, it does not quantitatively match the shape 
of measured speed vs d c . Below d c « 0.5, the model pre- 
dicts that v x (d c ) increases monotonically with increasing 
d c . This trend is in accord with the experimental ob- 
servations; however the model prediction is consistently 
higher than measured speed. Above d c « 0.5, the model 
predicts that v x (d c ) is independent of d Cl but the mea- 
sured speed decreases with increasing d c and is lower than 
the model prediction. We now discuss possible reasons 
for the discrepancies at low d c < 0.1 (labeled region I in 
Fig. 6b), intermediate 0.1 < d c < 0.5 (region II) and high 
d c > 0.5 (region III). 

In region I, d c is small so that uj s and thus a become 
large enough to ensure that the robot is in the swimming 
mode (i.e. movement without solidification). Here the 
model assumes rotary walking and predicts zero speed. 
In experiment, the robot can still advance slowly at each 
step due to thrust forces from continuously slipping limbs 
generated by frictional drag |25j and/or inertial move- 
ment of material. This thrust competes with friction 
from belly drag, and as in [TT], we find that these re- 
sult in low average speed of v x « 1 cm/sec. 

The model's overestimate of speed at high d c (region 
III) is a result of tripod overlap and can be readily under- 



stood. When there is no tripod overlap (only one tripod 
with ground contact at any given time) each tripod ad- 
vances the robot a distance s for a total displacement of 
2s per period. This is the case for d c < 0.5. However in 
the limit of d c — 1 both tripods are simultaneously in the 
slow phase as the duration of the fast phase is zero. The 
simultaneous slow phases generate a total displacement 
of just s [26] instead of 2s without tripod overlap. As 
a result, the predicted speed at d c = 1 must be halved 
(i.e. v x = J; = 1^). Lacking a way to quantify the tri- 
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FIG. 6. Comparison of data and anisotropic force model pre- 
dictions for v x (8 s ) &ndv x (d c ). (a) Measured v x (9 S , 8o — —0.5) 
(blue circles) through the speed maximum (Fig. 2a) deviates 
from the anisotropic force model prediction (blue curve) at 
large 9 S due to limb overlap, (b) Measured v x (d c ) (blue cir- 
cles) for {9 a = 1.5,0o = —0.5} is maximum at d c ~ 0.5. 
The model (solid blue curve) accurately predicts the speed 
for d c ~ 0.5 but is inaccurate elsewhere due to contributions 
of swimming neglected by the model (region I), a decrease in 
rotary walking range from cratering induced depth reduction 
and unequal penetration forces developed by c-legs on oppo- 
site sides of the body (II), and tripod overlap at high d c (III). 
Model prediction with tripod overlap included (dashed blue 
line) better matches the data in region III. 



pod overlap effect, we assume that the reduction of step 
length from 2s to s is linear with d c for d c > 0.5; the 
data is in good agreement with this prediction (dashed 
blue curve in Fig. 5) . This reduction in speed is a purely 
kinematic effect that is inherent to the rotary walking 
gait at high d c . 

Two plausible mechanisms explain the model's over- 
estimate of speed at intermediate d c (region II): hole 
digging and uneven weight distribution. Lateral obser- 
vations of the robot kinematics at low d c show that the 
rapid motion of the c-leg during the slow phase throws 
significant numbers of particles out of the limb's path 
which creates a depression. For a deep enough hole, ro- 
tary walking is impossible due to the reduced penetration 
depth of the leg below the now lower surface of the de- 
pression. The second mechanism concerns the model's 
assumption of uniform weight distribution between the 
three legs of the tripod. In the d c range where the robot 
advances slower than predicted, observations show that 
the robot rotates in the horizontal plane. Rotation oc- 
curs in this transition region between pure swimming and 
pure rotary walking because the side of the robot with 
two c-legs in ground contact undergoes rotary walking 
while the opposite side is in the swimming mode. Due 
to the increased gravitational and inertial forces on the 
single c-leg, the penetration forces are never sufficient to 
achieve ground solidification under the leg even at the 
maximum penetration depth. 

CONCLUSIONS 

We have built upon our previous experiments and 
models of a legged robot, SandBot, to explore how 
changes in limb kinematics affect locomotion on granular 
media. We found that even when moving on controlled 
granular media of fixed volume fraction at fixed cycle- 
averaged limb frequency, speed remains sensitive to vari- 
ations in gait parameters that control angular extent, an- 
gular location, and temporal duty factor of the slow phase 
of the limb cycle. We showed that the assumptions in a 
previously introduced model (which accurately predicted 
speed as a function of limb frequency and volume frac- 
tion) had to be modified to incorporate an anisotropic 
penetration force during rotational intrusion into gran- 
ular media as well as changes in acceleration of the leg 
as gait parameters were varied. With these modifica- 
tions the model was able to capture speed as a function 
of angular extent and angular location. The model also 
indicates that as duty cycle is changed, effects due to si- 
multaneous limb pairs (tripods) in ground contact, rapid 
limb impact into sand, and unequal weight distribution 
on limbs become important. 

Our experiments and modified model explain why gait 
parameters that allow the robot to rapidly bounce over 
hard ground lead to loss of performance on granular me- 



dia. They demonstrate how the angular extent and lo- 
cation of the slow phase must be adjusted to optimize 
interaction with granular media by minimizing inertial 
force and limb interference, and maximizing the use of 
solid properties of granular media. Further studies of 
SandBot guided by our kinematic model should reveal 
how physical parameters of both robot (mass distribu- 
tion, limb compliance, limb shape, belly shape) and the 
environment (grain friction, density, incline angle, grav- 
ity) control the solid-fluid transition and thus affect the 
limb-ground interaction and performance. However, ad- 
vances are required in theory and experimental charac- 
terization of complex media. Otherwise we must continue 
to rely on empirical force laws specific to particular ge- 
ometries, kinematics and granular media. 

The existence of a speed optimum in gait parameter 
space implies that control of limb kinematics is critical 
to move effectively on granular media, whether actively 
through sensory feedback, and/or passively through me- 
chanical feedback. Future work should compare these 
results to investigations of gait optimization on hard 
ground |13j . The differences in limb kinematics on sand 
compared to hard ground are intriguing because on hard 
ground performance is optimized by making the robot 
bounce. However, this carries with it the risk of yaw, 
pitch and roll instability due to mismanaged kinetic en- 
ergy. On granular media such instabilities appear rare; 
instead most gait parameters (see Fig. 3a) result in little 
or no forward movement due to mismanaged fluidization 
of the ground. Thus, our results could have a practical 
benefit as they suggest strategies for improving the per- 
formance of current machines [27- 29 on variable terrain 
via new limb and foot designs and control strategies. 

Finally, an enormous number of organisms contend 
with sand [30| . moving on the surface (or even swim- 
ming within it [31]). While the observed phenomena 
and proposed locomotion modes (e.g. rotary walking) 
appear specific to SandBot and its c-shaped limbs, the 
underlying principles could apply to locomotion of or- 
ganisms on yielding substrates. For example, our re- 
cent work on terrestrial hatchling sea turtle locomotion 
demonstrates that their effective movement on sand pro- 
ceeds through solidification of the granular medium [32) . 
Integrated studies of biological organisms and physical 
models can provide hypotheses [S3] for passive and ac- 
tive neuro-mechanical [34] control strategies as well as 
better understanding of energetics |35j for movement on 
complex terrain. 
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